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I  INTRODUCTION 


w 


This  report  is  the  thirteenth  bimonthly  progress  report  covering 
SRI  International’s  current  study  of  the  response  of  adaption  kits  in 
accidental  side  impacts.  The  primary  goal  of  the  project  is  the  further 
development  of  finite  element  models  for  a  generic  missile  structure  so 
that  adaption  kit  response,  particularly  component  accelerations,  can  be 
predicted  for  side  impacts  of  the  missile.  An  additional  goal  is  experi¬ 
mental  determination  of  the  impact  response  of  a  missile  structure  that 
contains  a  hard  link  safe-arm  device. 


The  work  planned  for  this  contract  consists  of  five  tasks.  In 
Task  1,  we  are  further  developing  finite  element  models  for  substructures 
and  checking  the  models'  predictions  against  experimental  data.  In 
Task  2,  we  are  using  the  models  to  calculate  the  response  of  complete 
structures,  again  by  comparing  the  predictions  with  experimental  results 
insofar  as  feasible.  In  Task  3,  we  are  transferring  the  analytical  capa¬ 
bilities  developed  in  Tasks  1  and  2  by  implementing  the  SUPER  code  at 
AARADCOM.  In  Task  4,  we  are  performing  impact  experiments  on  scale  models 
of  a  specific  prototype  missile  structure  to  help  plan  some  full-scale 
tests.  In  task  5,  we  are  developing  a  new  bolt  element  for  use  in  three- 
dimensional  calculations  ./^For  a  more  detailed  description  of  Tasks  1 
through  3,  refer  to  Bimonthly  Progress  Report  No.  1,  July  1978.  Task  4 
is  described  fully  in  Bimonthly  Progress  Report  No.  3,  November  1978. 

Task  5  has  only  recently  been  added  and  is  described  more  fully  in  the 
next  paragraph. 

Task  5  is  to  develop  a  new  bolt  element  for  use  in  three-dimensional 
calculations.  The  new  bolt  element  will  allow  the  SUPER  finite  element 
code  to  predict  the  response  of  adaption  kits  in  side  impacts  in  which 
plate-bending  occurs.  The  bolt  element  developed  previously  can  be  used 
only  for  calculations  of  in-plane  response  because  it  transmits  only  in¬ 
plane  forces.  The  new  bolt  element  will  use  a  nonlinear  spring  with 
modifications  for  unloading  and  for  breaking.  The  bolt  element  formula¬ 
tion  will  be  included  in  the  final  verson  of  the  SUPER  code,  which  will 
be  delivered  to  AARADCOM. 
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II  PROGRESS 


The  new  bolt  element  (Task  5)  was  developed  during  this  reporting 
period . 

Bolt  Element  Configuration 

The  purpose  of  the  bolt  elements  used  in  the  SUPER  calculations  is 
to  simulate  bolt  shearing  and  bolt  hole  stretching  with  a  single  element. 
Standard  elements  are  used,  but  the  material  properties  of  these  elements 
are  chosen  so  that  the  shear  force  versus  shear  displacement  relation  of 
the  bolt  element  approximates  the  results  of  static  lap  shear  tests. 

The  new  bolt  element  is  modeled  with  a  nonlinear  spring  that  trans¬ 
mits  only  axial  force.  Figure  1  shows  the  layout  of  the  bolt  element 
for  modeling  bolts  that  attach  the  AK  plate  to  the  ring.  (The  bolt 
element  can  also  be  used  to  model  bolts  that  attach  AK  components  to 
the  AK  plate.)  The  horizontal  component  of  the  axial  force  in  the 
spring  element  represents  the  in-plane  shear  force  that  is  transmitted 
between  the  plate  and  ring  (or  between  the  plate  and  the  AK  components). 
The  relative  motion  between  the  plate  and  the  ring  due  to  impact  will 
produce  a  relative  horizontal  displacement  d  at  the  ends  of  the  spring 
element  as  shown  in  Figure  2(a).  This  relative  displacement  results  in 
an  axial  stretch  in  the  spring  element,  as  shown  in  Figure  2(b).  Asso¬ 
ciated  with  this  axial  stretch  is  an  axial  force  that  has  a  horizontal 
component.  The  relationship  between  the  axial  force  and  axial  stretch 
is  chosen  so  that  the  horizontal  force  versus  horizontal  displacement 
approximates  the  results  of  the  lap  shear  tests. 

SUPER  Nonlinear  Spring 

The  finite  element  used  for  the  bolt  is  a  nonlinear  spring.  Figure  3 
shows  the  axial  force  N  •  £(A.i)  versus  change  in  length  M.  for  a  general 
nonlinear  spring.  In  the  SUPER  code,  this  curve  is  approximated  by  three 
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FIGURE  2  RELATIVE  DISPLACEMENT  PROOUCED  AT  THE  ENDS  OF  THE  SPRING 
BY  IMPACT  MOTION 
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straight  line  segments  as  in  Figure  4.  The  slopes  of  the  straight  lines 
are  the  stiffnesses  of  the  spring.  The  total  axial  force  in  the  spring 
between  AJl^  and  A£^+^  is  given  by 

N  •  [Z  Ki(A£i  -  A£i_1)  ]  +  Ki+1(A£i+1  -  AJ^) 

Figure  5  shows  the  force  and  displacement  of  the  bolt  element.  In 
the  figure,  Z^  is  the  initial  length  of  the  element,  d  is  the  relative 
horizontal  displacement  of  the  ends  of  the  element,  N  is  the  axial  force, 
F  is  the  horizontal  component  of  the  axial  force,  and  9  is  the  angle  of 
rotation  of  the  element  from  its  original  position.  The  change  in  length 
A Z  is  given  by 

A Z  -  y/d2  +  l2  -  Z 

o  o 

The  axial  force  is  a  function  of  AJl  as  indicated  by 
N  -  f(A Z) 

The  horizontal  component  of  N  is 

F  *  N  sin  0 

-  f  (AJl)  sin  9 

where 


sin  0  -  . “ . - 

y d2  +  zz 

o 

As  indicated  in  Figure  4,  the  spring  element  in  the  SUPER  code  can  have 
three  stiffnesses.  Thus,  equations  for  horizontal  force  F  versus  hori¬ 
zontal  relative  displacement  d  are,  in  general: 
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FIGURE  4  STRAIGHT-LINE  APPROXIMATION  OF 
AXIAL  FORCE-OISPLACEMENT  FOR 
A  NONLINEAR  SPRING 
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FIGURE  5  AXIAL  FORCE-DISPLACEMENT 
OIAGRAM  FOR  BOLT  ELEMENT 
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K,(v/d2  +  iz  -  l  )  sin  0 
o  o 


for 


(KiA^  +  K2(\/d2  +  £2  -l  -  AJlj)]  sin  6  for 


+  K2(A£2  -  Adi)  +  K3('/d2  +  £2  -  l  - 


o  o 

for 


Ail  <  A£i 

A£j  <_  Ail  <  A£2 

A£2) ]  sin  9 
A£2  £  Ail  <  Ail  3 


Loading  Characteristics 

The  lap  shear  test  results  suggest  that  the  horizontal  component  of 
the  axial  force  versus  horizontal  displacement  resembles  a  sine  function 
of  the  form 

F  *  A  sin  <j> 

where  A  is  an  arbitrary  constant  and  $  is  a  parameter  in  the  range 
o  <  $  <  tt/2.  Thus,  a  good  approximation  to  the  test  results  may  be 
obtained  by  using  a  spring  with  a  constant  axial  force  Nj,  so  that 

F  *  Nj  sin  9 

Since  the  spring  element  is  governed  by  stiffnesses  rather  than  by  forces, 
we  use  a  spring  with  an  initial  stiffness  Kj  up  to  a  constant  force  Nj 
and  K2,  K3  equal  to  zero.  Choosing  *  790,866  N/mm,  iQ  *  0.381  mm, 

A£i  *  0.00132  mm  gives  a  good  fit  to  the  lap  shear  data. 

Referring  back  to  the  equations  for  the  horizontal  force,  we  have 

F  »  Ki  (\/d2  +  £2  -  £  )  sin  9  for  Ail  <  0.00132  mm 

1  o  o 

F  ■  Ki  A£i  sin  9  for  Ail  >_  0.00132  ram 


and 


-  Ki  AS.!  -  1044  N 
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Figure  6  shows  che  axial  force-change  in  length  relationship  for  the 
spring  element.  A  plot  of  the  shear  force  versus  shear  displacement  for 
the  bolt  element  model  and  the  lap  shear  test  results  is  shown  in  Figure  7. 

Unloading  Characteristics 

Figure  8  shows  the  axial  force,  its  horizontal  component,  and  the 
relative  displacement  of  the  bolt  element  under  unloading;  di  and  d2  are 
any  two  points  along  the  unloading  curve  of  the  bolt  element  and  di  >  d2 . 

The  change  in  length  of  the  spring  element  is  given  by 

M  -  \/di2  +  i  2  -  \/do2  +  l  2 
u  o  c  o 

Provided  the  spring  is  not  loaded  beyond  the  rupture  point,  the  change 
in  axial  force  AN  associated  with  the  change  in  length  A£  is  assumed 
to  be  linear,  given  by 


AN  *  K  Ax. 

u  u 

where  is  the  coefficient  of  stiffness  for  unloading  of  the  spring 
element.  The  horizontal  component  of  the  axial  force  is 

F  -  (N  -  AN)  sin  9 
di 

F;  *  N - - 

\/d  i 2  +  x  2 
1  o 


d2 


A  spring  element  with  *  8584  N/mm  will  give  us  an  unloading  curve  for 
the  horizontal  component  of  the  axial  force  versus  relative  horizontal 
displacement  that  matches  the  experimental  results.  A  comparison  of  the 
bolt  element  model  with  the  lap  shear  test  result  for  both  loading  and 
unloading  was  shown  in  Figure  7. 
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FIGURE  6 
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AXIAL  FORCE-CHANGE  IN  LENGTH  RELATION¬ 
SHIP  FOR  THE  SPRING  ELEMENT 
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FIGURE  7  BOLT  ELEMENT  MODELED  AS  A  NONLINEAR  SPRING 


FIGURE  8 


FORCE-DISPLACEMENT  RELATIONSHIP 
OF  THE  BOLT  ELEMENT  UNDER 
UNLOADING 
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Ill  FUTURE  WORK 


During  the  next  reporting  period,  a  SUPER  calculation  of  a  three- 
dimensional  multicomponent  structure  will  be  conducted.  The  calculation 
will  test  the  adequacy  of  the  new  bolt  element  for  use  in  calculations 
involving  plate  bending.  The  results  of  the  calculation  will  also  be 
compared  qualitatively  with  an  experiment  on  a  complete  missile  section. 

( 

IV  FINANCIAL  STATUS 

As  of  27  June  1980,  $84,500  has  been  spent  on  labor  (1671  supervisory 
and  professional  hours,  968  technical  and  clerical  hours),  and  $11,631 
has  been  spent  on  materials  and  services.  Of  the  total  contract  funds 
($98,435),  the  balance  remaining  is  $2,304.  A  performance  and  cost  report 
for  this  period  is  attached. 


PERFORMANCE  AND  COST  REPORT 
DAAK-10-78-C-0158 
PYU-7422  -  Report  No.  13 


Reporting  Period:  3  May  1980  -  2  July  1980 

Hours 

Total  hours  expended  to  date: 

Supervisory  and  Sr.  Professional  Personnel 
Professional  Personnel 
Technician  and  Clerical 

Cumulative  total  hours  to  date 
Percent  of  total  hours  expended  to  date 

Funds 

Funds  expended  during  the  reporting  period 

Funds  expended  to  dace 

Percent  of  total  funds  expended  to  date 

Work 

Percent  of  work  completed  during  the  reporting 
period 

Percent  of  work  completed  to  date 
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275 

1396 

968 

2639 

98% 


$  2,235 
$96,130 
98% 


3% 

98% 
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